Quantum control and entanglement using periodic driving fields 
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We propose a scheme for producing directed motion in a lattice system by applying a periodic 
driving potential. By controlling the dynamics by means of the effect known as coherent destruction 
of tunneling, we demonstrate a novel ratchet-like effect that enables particles to be coherently 
manipulated and steered without requiring local control. Entanglement between particles can also 
be controllably generated, which points to the attractive possibility of using this technique for 
quantum information processing. 
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Introduction Controlling the time-evolution of quan- 
tum states, and engineering entanglement between quan- 
tum particles, are two of the major tasks required for 
quantum information processing. Although the majority 
of experimental demonstrations of entanglement distri- 
bution have so far employed photons, future practical 
implementations of quantum computers will almost cer- 
tainly be based on condensed matter systems. Efforts 
in this direction include controlling electronic charge or 
spin ( "spintronic" ) degrees of freedom in coupled quan- 
tum dots [l|, manipulating the dynamics of Joscphson 
junctions [2J, and using spin chains [3| as quantum com- 
munication channels. Recently bosons held in optical 
lattice potentials have also been suggested as possible 
candidates, using controlled collisions Q or the dipole- 
dipole interaction [B[ to manipulate the system. 

In this work we demonstrate how particles in a lat- 
tice potential can be controlled by applying an oscilla- 
tory driving field. Control is achieved by using the effect 
termed "coherent destruction of tunneling" (CDT) Q 
to generate a ratchet-like motion. Unlike the majority 
of ratchets which depend on the presence of dissipation 
to operate, this motion results solely from the breaking 
of spatial and time symmetries in the system, and so 
preserves the quantum coherence of the system. Using 
this scheme, a pair of particles can be brought together 
and allowed to interact - thereby becoming entangled 
- and then subsequently separated again to create en- 
tanglement between distinct spatial locations. We are 
thus able to both selectively move and entangle quantum 
particles using only the global operation of varying the 
parameters of the driving field, thus avoiding any need 
for the individual addressing of lattice sites. 

Model The specific physical system that we consider 
consists of ultracold bosonic atoms, confined in a one- 
dimensional optical lattice potential created by the su- 
perposition of counter-propagating laser beams. This 
form of confinement provides an extremely clean and con- 
trollable lattice potential, and in addition their high de- 
gree of isolation from the environment gives these sys- 
tems rather long decoherence times, making them ideal 



for studying quantum coherent phenomena. 

The single-particle dynamics of the system can be de- 
scribed extremely well by the Hamiltonian 0] 



H = 



Jij a\a j 



H.c. 



K sin cjt 



(1) 



where Ou i ja\ are the standard bosonic destruc- 
tion/creation operators and the tunneling amplitudes 
Jij connect nearest-neighbor sites Without loss 



of generality we shall henceforth take J, : 



J. The 



amplitude and frequency of the time-dependent driving 
field are described by the parameters K and to, and 
Xi is the spatial location of the ith lattice site. This 
form of linear potential has already been used in cold 
atom experiments [8] to induce CDT, and can be 
straightforwardly implemented in an optical lattice by 
introducing a periodic phase-modulation to one of the 
laser fields providing the standing wave potential. 

Since the Hamiltonian |T]) is invariant under discrete 
translations in time of the drive-period, H(t) = H(t + 
nT), the Floquet theorem allows us to write solutions of 
the Schrodinger equation as \<f> n {t)) = \u n (t)} exp [— ite n ], 
where \u n (t)) is a T-periodic function called the Floquet 
function, and e n is termed the quasienergy. When two 
quasienergies approach degeneracy, the timescale to ob- 
serve tunneling between the associated Floquet states di- 
verges, and accordingly the tunneling between them ap- 
pears suppressed. In the limit of high-frequency (when u) 
is the dominant e nerg y-scale of the problem) it may be 
indeed shown 0, 0, 1 1 OL fxxj] that the driven system behaves 
like the undriven one, but with renormalized tunneling 
amplitudes. For sinusoidal driving this renormalization 
takes the form J e g = J^Jo(Kx/uj), where j7o is the zeroth 
Bessel function of the first kind and x is the intersite sep- 
aration. Thus when Kx/lu is equal to a zero of j7o the 
system's tunneling dynamics are frozen, producing CDT. 

This way of regulating the tunneling between sites 
has been recently proposed to control the Mott-insulator 
transition in Bosc-Einstcin condensates 12, la. To ob- 



tain a ratchet effect, however, it is necessary to distin- 
guish between motion to the left and motion to the right. 
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FIG. 1: We consider a bipartite lattice with two spacings: xi 
and X2- Permitted tunneling processes are shown by the ar- 
rows, (a) In the absence of a driving field, a particle initialized 
in a given lattice site (the filled circle) can tunnel to either of 
its neighbors (empty circles). From there it can spread over 
the entire lattice, (b) If the lattice is sinusoidally-driven with 
a field such that Jo(Kxi/lj) = then the tunneling processes 
between sites separated by x\ are suppressed, and the particle 
can only tunnel to its right neighbor, (c) Conversely, if the 
driving field satisfies Jo[KxiliJ) = 0, then tunneling is only 
permitted between sites separated by xi, and the particle can 
only tunnel to its left neighbor. 



This may be achieved by noting that the argument of 
the Bessel function depends on the potential difference 
between neighboring sites, and thus on their spatial sep- 
aration. Accordingly we consider a bipartite lattice of 
form ABABAB, as shown in Fig[TJ in which the AB sep- 
aration is not equal to that between BA. One possible 
realization of this would be a chain of coupled double- 
well potentials [3] • Consider placing a single particle in 
the center of this lattice (FigOJi). In the absence of a 
driving field it will rapidly disperse by tunneling to both 
its neighbors. If, however, the system is driven by a high- 
frequency sinusoidal potential such that Jo(Kxi/ui) = 0, 
then tunneling between sites separated by x\ is destroyed 
and the lattice divides into a set of disconnected dimers 
(AB)(AB)(AB). In this case the particle is unable to 
spread over the lattice, and is restricted to making a 
Rabi oscillation (FigQJ>) between its initial location and 
its neighbor to the right. The frequency of this oscillation 
is determined by the value of the renormalized tunneling 
for this process J e ff = JJo(Kx2/uj), which will in general 
be non-zero. Conversely, if the parameters of the driving 
field are chosen such that Jo {Kxz/uj) — then the lattice 
dimerises as {BA)(BA)(BA) as shown in Fig[TJ;, and the 
Rabi oscillation will occur between the initial site and its 
neighbor to the left. 



Results To verify this effect, we show in Fig [2^ the re- 
sults of a numerical simulation of a single particle in a 
16-site system with x\ — \ and x 2 — 0.75. The frequency 
of the driving is set to a high value of u> = 32 J to ensure 
that the system is in the high-frequency regime, while its 
amplitude satisfies Kxi/cu = 2.4048 - the first zero of 
J7b- As the analysis predicts, the particle indeed simply 
oscillates between its initial location and one of its neigh- 
bors (that separated by a distance of X2), since tunneling 
between sites separated by x\ has been suppressed. 

The crispness of the Rabi oscillation immediately sug- 
gests a scheme to produce directed motion. If we denote 
the period of this oscillation by T2, then at t = T2/2 the 
particle has completely tunneled from its initial location 
i to its neighbor i + 1. If at this time the parameters of 
the field are altered so that Kx-ijoj = 2.4048, then this 
tunneling process is suppressed and the particle instead 
begins to make a Rabi oscillation with period T± between 
sites i + 1 and i + 2. When a time interval of Ti/2 has 
elapsed the particle has completely tunneled to site i + 2. 
The driving field can then be switched back to its original 
values and the procedure repeated. This has the effect of 
stepping the particle through the lattice in a sequence of 
discrete moves. 



An example of a driving field that can produce this 
effect is shown schematically in Fig[3] It can be thought 
of as a high-frequency "carrier wave" whose amplitude 
is modulated by a squarewave envelope. The lower am- 
plitude segment suppresses tunneling between sites sep- 
arated by xi and has a duration such that the particle 
tunnels exactly to its other neighbor (separated by £2): 
the reverse is true for the higher amplitude segments. 
FigEJa shows the response of the single-particle system 
to this field. Instead of the two-site Rabi oscillation seen 
previously, the particle now advances to the left in a se- 
ries of well-defined steps. Conversely, if the order of the 
modulation is interchanged, the particle will propagate 
solely to the right as shown in Fig[2j;. It is interesting 
to note that the direction of the particle's propagation 
also depends on which site of the double-well it is ini- 
tialized; particles started in the left well will move in the 
opposite direction to those placed in the right. The di- 
rection of motion thus depends on both the parity of the 
lattice site and the order of modulation, in a way not 
seen in standard dissipative ratchets. This flexibility re- 
quires, however, excellent control over the localization of 
the initial state. 

We have so far considered the extreme cases in which 
propagation only occurs in one direction. If, however, 
the duration of the initial modulation is not exactly half 
a Rabi-period, the initial tunneling process will not be 
complete. Consequently the particle will divide into two 
parts, and under the subsequent influence of the driving 
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FIG. 2: Simulation of a single boson in a 16-site system, (a) 
Under a periodic driving field of amplitude Kxi/u) = 2.4048 
we obtain the situation illustrated in Fig[lp - tunneling be- 
tween sites separated by x\ is suppressed. Consequently the 
particle makes a Rabi oscillation to just one of its neighbors, 
(b) Under the driving field shown in Fig[3] the xi and xi 
tunneling processes are periodically opened and closed, pro- 
ducing a ratchet-like motion of the particle, (c) Interchanging 
the order of the modulation of the driving field produces mo- 
tion in the opposite direction, (d) Shortening the duration of 
the initial modulation means that the first tunneling process 
will be incomplete. If the initial modulation time is halved 
the particle splits into two equal parts, and under the driving 
field each part propagates in different directions. 

field one part will propagate to the left while the other 
moves to the right. In FigJ^Ji we show that if the ini- 
tial modulation has a duration of Ti /4 the particle splits 
in half. In this way, the driving field can not only be 
used to control the motion of a particle, but also as a 
quantum beam splitter to divide a particle into a given 
superposition of left and right propagating components. 



We have so far just considered single-particle manipu- 
lation. However, the ability to control the inter-site tun- 
neling also enables us to entangle particles and thereby 
realize quantum gates To demonstrate this, we first 
consider a two-site model occupied by two particles (o 
and b) that are distinguishable. This can be realized, for 
example, by using bosonic atoms 0, 0] with two different 
internal states [15| . For simplicity, we model the interac- 
tion between the particles as a Hubbard interaction 

A> = -1 ) > ( 2 ) 

i 

where U sets the interaction-strength and m = nf + rq 
gives the total number of bosons occupying site i. The 



FIG. 3: The driving field producing the ratchet-like motion 
seen in Fig[2]consists of a high-frequency sinusoidal oscillation 
with a modulated amplitude. The two values of the modula- 
tion satisfy K/ui — 2.4048/a;i, and act to suppress the tun- 
neling between sites separated by distances of Xi respectively. 
The time-intervals T\ and Tb are the Rabi periods for the 
renormalized tunneling processes between sites separated by 
xi and x%. 



dynamics of this system is governed by the interplay be- 
tween the kinetic energy and the interaction, and conse- 
quently can exhibit a rather complicated time-evolution. 
If, however, U is much larger than the tunneling ampli- 
tude, the ground state of the system will then approx- 
imately consist of each site holding one particle, from 
which the doubly-occupied states will be separated by 
an energy gap of ~ U. In this case canonical perturba- 
tion theory can be applied to eliminate the higher en- 
ergy states, with the result that the Hubbard interaction 
maps to an effective Heisenberg term, with exchange con- 
stant given by Jjr — ^Jgjj/?/. This mapping considerably 
simplifies analysis of the system's dynamics, and reveals 
that if the system is initialized in the state \a, 6), the two 
particles swap positions after a time-interval t — it/ Jjj, 
while after ts = t/2Jh the maximally entangled state, 
(|a, b) + i\b,a))/y2, is produced. Applying the Heisen- 
berg interaction for a duration of ts thus realizes the 
VSWAP operation. 

In Fig[4]we show the time-evolution of a 10-site lattice, 
initialized with a boson of type a in the first site, and a 
boson of type b in the last. To ensure the validity of the 
mapping to the Heisenberg interaction, we require a high 
value for U. However, the time required to entangle the 
particles is proportional to U, and so to complete as many 
quantum gate operations as possible within the system's 
decoherence time we would like to take U to be as small 
as possible. We thus consider an intermediate value of 
U = 4J, and as before use a driving frequency of lu = 
32J to place the system in the high-frequency regime. 
Under the influence of the driving field the two bosons 
are progressively stepped through the lattice towards the 
two central sites, whereupon the amplitude of the driving 
is held at a constant value to retain the two particles 
there. After being held there for a time interval of tg, 
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FIG. 4: Above: Time development of the system initialized 
with one boson in the first lattice site, and another in the final 
site. Under the influence of the driving field they are moved 
toward the center of the array, held there for a time-interval 
t$, and then restored to their original positions. Below: The 
pairwise concurrence present in the system. Initially this is 
zero since the initial state is not entangled, but rises as the 
two bosons interact (via Heisenberg exchange). The interval 
over which the interaction occurs is chosen to maximize the 
concurrence, by producing a maximally entangled state. 

the particles are then separated and returned to the first 
and last lattice sites. 



From FigHk it can be clearly seen that the parti- 
cles remain highly localized in space, and accordingly 
the probability distribution is generally peaked at two 
sites. This permits a simple measurement of the entan- 
glement present in the system by projecting the wave- 
function onto just these two sites, and then evaluating 
the two-particle concurrence C(t) [lj]], plotted in FigHb. 
Initially the concurrence is zero since the particles have 
not interacted, and so the two-particle wavefunction is 
factorisable. This remains true as the particles approach 
each other, until they reach the two central sites. Driven 
by the Heisenberg interaction, the concurrence then rises 
from zero following the approximate time-dependence 
C(t) = | sin Jgt\. The high-frequency ripples visible in 
this quantity arise from the influence of the higher energy 
states: if U is increased these ripples will be quenched, 
but equally Jjr will be reduced, and so the time-scale 
for entanglement to occur will increase. When the par- 
ticles are separated the degree of entanglement remains 



"frozen-in" at its final value. It is thus possible to gen- 
erate any desired degree of entanglement by controlling 
the period during which the particles interact. When this 
period is equal to tg, as shown in FigHl the entanglement 
is maximized, and the final state of the system thus rep- 
resents a mesoscopically separated, maximally-entangled 
two-particle state. 

Conclusions In summary, we have shown how a peri- 
odic driving field can induce a novel ratchet-like motion, 
which can be employed to selectively guide and divide 
particles. In addition the interaction between the parti- 
cles can be used to entangle them, and thus realize fun- 
damental two-qubit quantum gates, such as -^/SWAP . 
In this work we have just considered a one-dimensional 
geometry, but an attractive aspect of optical lattices is 
the possibility of generating a higher-dimensional flij lat- 
tice potentials, which would allow parallel processing of 
qubits, and thus be used to greatly enhance their error 
tolerance. Finally although we have specifically consid- 
ered a system of ultracold bosons, the method we have 
described could equally be applied to optically-confined 
fermionic atoms [18j |. or to electronic transport in sys- 
tems such as coupled quantum dots (l9j or molecular 
wires [H. 
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